*** Which type of branches ****
use "$Data_generated/Commune.dta", replace

egen sum_meetings=sum(AF_Meetings) if year>1924 & year<1929, by(year) 

egen sum_sect=sum(nb_sect), by(year)
egen sum_sectpres=sum(nb_sectpres), by(year)
egen sum_sectperm=sum(nb_sectperm), by(year)
duplicates drop year, force
tsset year

gen max=1000
twoway (bar max year if inrange(year, 1926.5, 1927.5), color(gs12)) (line sum_sect sum_sectpres sum_sectperm year, xtitle("Year") ytitle("Nb of Branches") xlabel(1921(1)1936) lcolor(gs1 gs6 gs9) legend(order(2 "Nb of Branches" 3 "Nb of Branches with a President" 4 " Nb of Branches with Premises"))  scheme(stcolor_alt) text(600 1927 "Condemnation", orientation(vertical))) (scatter sum_sect  year if year>1920 & year<1936, mcolor(gs6) msymbol(T)) (scatter sum_sectpres year if year>1920 & year<1936, mcolor(gs6) msymbol(S)) (scatter sum_sectperm year if year>1920 & year<1936, mcolor(gs9) msymbol(X))
graph export "$Figures/Figure6a.pdf", name(Graph) replace
graph export "$Figures/Figure6a.png", name(Graph) replace



**** Figure 2: Evolution of number of Meetings over time ****
tsset year
replace max=10000
twoway (bar max year if inrange(year, 1926.5, 1927.5), color(gs12)) (line sum_meetings year if year>1924 & year<1929, xtitle("Year") ytitle("Nb of Branches") xlabel(1925(1)1928) lcolor(gs1 gs6 gs9) legend(off) text(8000 1927 "Condemnation", orientation(vertical))) (scatter sum_meetings  year if year>1924 & year<1929, mcolor(gs6) msymbol(T))
graph export "$Figures/Figure6b.pdf", name(Graph) replace
graph export "$Figures/Figure6b.png", name(Graph) replace
